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Abstract. We study b, c quark hadronization from QGP. We obtain the yields of charm and bottom flavored 
hadrons within the statistical hadronization model. The important novel feature of this study is that we 
take into account the high strangeness and entropy content of QGP, conserving strangeness and entropy 
yields at hadronization. 



PACS. 25.75.Nq Quark deconfmement - 12.38. Mh 
25.75.-q Relativistic heavy-ion collisions - 24. 10. Pa 



1 Introduction 

A relatively large number of hadrons containing charmed 
and bottom quarks are expected to be produced in heavy 
ion (AA) collisions at the Large Hadrons Collider (LHC). 
Because of their large mass c, c, b, b quarks are produced 
predominantly in primary parton-parton collisions pQ, at 
RHIC [2], and thus even more so at LHC. These heavy fla- 
vor quarks participate in the evolution of the dense QCD 
matter from the beginning. In view of the recent RHIC 
results it can be hoped that their momentum distribution 
could reach approximate thermalization within the dense 
QGP phase [3]. 

In our approach we will tacitly assume that the follow- 
ing evolution stages are present in heavy-ion collisions: 

1. Primary partons collide producing c, b quarks; 

2. A thermalized parton state within r = r t h — 0.25 — 
lfm/c is formed. By the end of this stage nearly all 
entropy is produced. 

3. The subsequent chemical equilibration: diverse ther- 
mal particle production reactions occur, allowing first 
the approach to chemical equilibrium by gluons g and 
light non-strange quarks q = u,d. 

4. The strangeness chemical equilibration within r ~ 5 
fm/c. 

5. The hadronization to final state near r ~ 10 fm/c. 

It is important to observe that in the presence of decon- 
fincd QGP phase heavy hadrons containing more than one 
heavy quark are made from heavy quarks created in differ- 
ent initial NN collisions. Therefore yields of these hadrons 
are expected to be enhanced as compared to yields seen 
in single NN collisions [4j[5] . We note that the Be (be, be) 
and J I ^ ice) and more generally all bound cc states yields 
were calculated before in the kinetic formation and disso- 
ciation models [IJ[6]. Our present work suggests that it is 
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important to account for the binding of heavy flavor with 
strangeness, an effect which depletes the eligible_supply 
of heavy flavor quarks which could form Bc(fec, be) and 
J/Hcc) 0. 

Enhanced production yield of multi-heavy hadrons can 
be considered to be an indicator of the presence of decon- 
fincd QGP phase for reasons which are analogue to those 
of multi-strange (anti) baryons 8J. Considering that we 
have little doubt that QGP is the state of matter formed 
in the very high energy AA interactions, the study of yields 
of multi-heavy hadrons is primarily explored in this work 
in order to falsify, or justify, features of the statistical 
hadronization model (SHM) employed or the model itself 
in the context of formation of the heavy flavor hadrons. 

For example, differing from other recent studies which 
assume that the hadron yields after hadronization are in 
chemical equilibrium 5,9 , we form the yields based on 
abundance of u, d, s quark pairs as these are available at 
the chemical freeze-out (particle formation) conditions in 
the quark-gluon phase. This approach is justified by the 
expectation that in a fast break-up of the QGP formed at 
RHIC and LHC the phase entropy and strangeness will be 
nearly conserved during the process of hadronization. We 
will investigate in quantitative terms how such chemical 
non-equilibrium yields, in the conditions we explore well 
above the chemical equilibrium abundance, influence the 
expected yields of single, and multi-heavy flavor hadrons. 

In the order to evaluate the yields of final state hadrons 
we enforce conservation of entropy, and the flavor s,c,b 
quark pair number during phase transition or transfor- 
mation. The faster the transition, the less likely it is that 
there is significant change in strange quark pair yield. Sim- 
ilarly, any entropy production is minimized when the en- 
tropy rich QGP breakup into the entropy poor HG occurs 
rapidly. The entropy conservation constraint fixes the fi- 
nal light quark yield. We assume a fast transition between 
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QGP and HG phases, such that all hadron yields are at 
the same physical conditions as in QGP breakup. 

In the evaluation of heavy particle yields we form ra- 
tios involving as normalizer the total heavy flavor yield, 
and for yields of particles with two heavy quarks we use 
as normalizer the product of total yields of corresponding 
heavy flavors such that the results we consider is as lit- 
tle as possible dependent on the unknown total yield of 
charm and bottom at RHIC and LHC The order of mag- 
nitude of the remaining dependence on heavy flavor yield 
is set by the ratio of yield of all particles with two heavy 
quarks to yield of particles with one heavy quark. This 
ratio depends on the density of heavy flavor at hadroniza- 
tion, (dN c /dy)/(dV/dy). The results we present for LHC 
are obtained for an assumed charm and bottom quark 
multiplicity: 



dNc 
dy 

dN b 
dy 



c = 10, 
6=1. 



(1) 
(2) 



and dV/dy = 800 fm 3 at T = 200 MeV. Theoretical 
cross sections of c and b quarks production for RHIC and 
LHC can be found in [TOUTTj . In certain situations we will 
explore how variation of the baseline yields Eq. ([TJ) and 
Eq. (|2|) impact the results. In particular among the yields 
of multi-heavy hadrons, this influence can be noticeable, 
see discussion in the end of section 16.11 We note that the 
number of b quarks can not change during expansion, be- 
cause of large mass >> T. It is nearly certain that all 
charm in QGP at RHIC is produced in the first parton 
collisions, for further discussion of LHC see Ref.[12] - it 
appears that for all practical purposes also in the more 
extreme thermal conditions at LHC charm is produced in 
the initial parton interactions. 

In order to form physical intuition about the prevailing 
conditions in the QGP phase at time of hadronization, we 
also evaluate the heavy quark chemical reference density, 
that is the magnitude of the chemical occupancy factor 
in QGP, considering the pre-established initial yields of c 
and b from parton collision. For this purpose we use in the 
deconfined QGP phase: 

m c = 1.2 GeV, 
m b = 4.2. GeV 

We also take Aj = l,i = u,d,s for all light flavors, since 
the deviation from particle-antiparticle yield symmetry is 
rather small and immaterial in the present discussion. 

When computing the yields of charmed (and bottom) 
mesons we will distinguish only strange and non-strange 
abundances, but not charged with non-charged (e.g. D~ (cd) 

with D (cu)). We assume that the experimental groups 
reporting results, depending on which types of D-meson 
were observed, can infer the total yield (charged+non- 
charged) which we present. We treat in similar way other 
heavy hadrons, always focusing on the heavy and the strange 
flavor content and not distinguishing the light flavor con- 
tent. 



Our paper is organized as follows: we first introduce 
the elements of the SHM model we use to evaluate heavy 
flavor hadron yields in section [2j This allows us to dis- 
cuss the relative yields of strange and non-strange heavy 
mesons in section [31 and we show how this result relates 
the value of the strangeness chemical (non-) equilibrium 
parameters. In this context, we also propose a multi-particle 
ratio as a measure of the hadronization temperature, and 
explore how a multi-temperature, staged, freeze-out would 
impact the relevant results. 

Before proceeding to obtain the main results of this 
work we introduce the notion of conservation of entropy in 
section[4]and strangeness in section[5l expected to be valid 
in the fast hadronization process at LHC, and discuss how 
this impacts the SHM statistical parameters. We consider 
the entropy in a system with evolving strangeness in sub- 
section 03] and show that the number of active degrees of 
freedom in a QGP is nearly constant. Another highlight 
is the discussion of sudden hadronization of strangeness 
and the associated values of hadron phase space parame- 
ters in subsection 15.41 Throughout this paper we will use 
explicitly and implicitly the properties of QGP fireball 
and hadron phase space regarding entropy and strangeness 
content developed in these two sections [4] and [5) 

We turn to discuss the heavy flavor hadron yields for 
given bulk QGP constraints in section [51 where we also 
compare when appropriate to the strangeness and light 
quarks chemical equilibrium results. We begin with a study 
of the charm and bottom quark phase space occupancy 
parameters (subsection [OJ and turn in subsection 16.21 to 
discussion of the yields of single heavy mesons, which we 
follow with discussion of yields of single heavy baryons in 
subsection 16.31 In last subsection 16.41 we present the ex- 
pected yields of the multi-heavy hadrons, in so far these 
can be considered in the grand canonical approach. We 
conclude our work with a brief summary in section [7l 



2 Statistical Hadronization Model (SHM) 

The statistical hadronization model arises from the Fermi 
multi-particle production model [13] . Fermi considered that 
all hadron production matrix elements are saturated to 
unity. This allows the use of the Fermi golden rule with 
the N-particle phase space to obtain the relative particle 
yields. In modern language this is SHM in micro-canonical 
ensemble, micro-canonical implies that discrete (flavor) 
quantum numbers and the energy are conserved exactly. 

The transition from micro-canonical to canonical, and 
grand-canonical ensembles simplifies the computational 
effort considerably [14] . This important step does not in 
our context introduce the hadron phase, although before 
the understanding of QGP this of course was the reaction 
picture: a highly compressed hadron gas matter evapo- 
rates particles. Today, it is the highly compressed hot 
quark-gluon matter that evaporates particles. In princi- 
ple, there is no necessity to introduce a hadron gas phase 
of matter in order to use SHM to describe particle pro- 
duction. 
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On the other hand, in order to understand the physical 
meaning of the parameters introduced to describe hadron 
phase space in grand-canonical ensemble, such as temper- 
ature T, it is quite convenient to imagine the existence 
of the hadron phase which follows the QGP phase. This 
can be taken to the extreme, and a long lasting, chemi- 
cally equilibrating phase of hadrons can be assumed, that 
follows in time the formation of the QGP fireball. Such 
a reaction picture may not agree with the fast evolving 
circumstance of a heavy ion collision. One should note 
that the study of hot hadron matter on the lattice within 
the realm of L-QCD involves at all times a fully equi- 
librated system. This will in key features differ from the 
non-equilibrium QGP properties accompanying the hadro- 
nization process. 

The important parameters of the SHM, which control 
the relative yields of particles, are the particle specific fu- 
gacity factor A and space occupancy factor 7. The fugacity 
is related to chemical potential /1 = Tin A. The occupancy 
7 is, nearly, the ratio of produced particles to the number 
of particle expected in chemical equilibrium. The actual 
momentum distribution is: 



d 6 N 
d 3 pd 3 x 



9 



1 



(2tt) 3 1 - 1 \- 1 e E / T ± 1 



7Ae 



-E/T 



(3) 

where the Boltzmann limit of the Fermi '(+)' and Bose 
'(— )' distributions is indicated, g is the degeneracy factor, 
T is the temperature and E = E(p) is the energy. 

The fugacity A is associated with a conserved quan- 
tum number, such as net-baryon number, net-strangeness, 
heavy flavor. Thus antiparticles have inverse value of A, 
and A evolution during the reaction process is related to 
the changes in densities due to dynamics such as expan- 
sion. 7 is the same for particles and antiparticles. Its value 
changes as a function of time even if the system does not 
expand, for it describes buildup of the particular particle 
species. For this reason 7 is changes rapidly during the 
reaction, while A is more constant. Thus it is 7 which car- 
ries the information about the time history of the reaction 
and the precise condition of particle production referred 
to as chemical freeze-out. 

The number of particles of type with mass rrii per 
unit of rapidity is in our approach given by: 



dNj _ cq dV 
dy 7jn * dy ' 



(4) 



Here dV/dy is system volume associated with the unit 
of rapidity, and n^ q is a Boltzmann particle density in 
chemical equilibrium: 



eq 

1 1 — 



9i l (|^3^ eX P(-V /p2+m ^ T) 



and 



Both, 



fp 
(2^) 

\i^9iW{mi/T), 



W{x) = x 2 K 2 {x) -> 2 for x -> 0. 



(5) 



(6) 



rriiC 



rrii, and kT — ► T, are measured in energy 



For the case of heavy flavors m » T, the dominant 
contribution to the Boltzmann integral Eq. ^ arises from 
p ~ V2mT, we do not probe the tails of the momentum 
distribution. Thus even when the momentum distribution 
is not well thermalized, the yield of heavy flavor hadrons 
can be described in therm of the thermal yields, Eq.Q, 
where: 



9 



T 3 



2tt 2 



2T 3 
15T 105 
8m" + 128 



exp(— m,/T) x 
2 



(7) 



Often one can use the first term alone for heavy flavor 
hadrons, since T/m << 1, however the asymptotic se- 
ries in Eq. ([7]) has limited validity. Our computations are 
all based on CERN recursive subroutine evaluation of the 
Bessel K 2 functions. 

We use occupancy factors 7^ and 7^ for QGP and 
hadronic gas phase respectively, tracking every quark fla- 
vor (i — q, s, b, c) . We assume that in the QGP phase 
the light quarks and gluons are adjusting fast to the ambi- 
ent conditions, and thus are in chemical equilibrium with 
7^ G — * 1. For heavy, and strange flavor, the value of 7^ at 
hadronization condition is given by the number of particles 
present, generated by prior kinetic processes, see Eq. (|4|). 

The yields of different quark flavors originate in differ- 
ent physical processes, such as production in initial col- 
lisions for c, b, s, and for s also production in thermal 
plasma processes. In general we thus cannot expect that 
^cb w ^ be near unity at hadronization. However, the ther- 
mal strangeness production process GG — > ss can nearly 
chemically equilibrate strangeness flavor in plasma formed 
at RHIC and/or LHC [12], and we will always consider, 
among other cases the limit 7^ — > 1 prior to hadroniza- 
tion. 

The yields of all hadrons after hadronization are also 
given by Eq. (U), which helps us to obtain 7? for hadrons. 
In general, the evaluation of hadron chemical parameters 
presents a more complicated case since they are composed 
from those of valance quarks in the hadron (three quarks, 
or a quark and anti quark). Therefore, in the coalescence 
picture, the phase space occupancy 7 H of hadrons will 
be the product of 7 H 's for each constituent quark. For 
example for the charmed meson D (cq) 



Id = 7 C % H 



(8) 



units when h,c,k—> 1 . 



To evaluate yields of final state hadrons we enforce 
conservation of entropy, and the flavor s,c, b quark pair 
number during phase transition or transformation. The 
faster the transition, the less likely is that there is signif- 
icant change in strange quark pair yield. Similarly, any 
entropy production is minimized when the entropy rich 
QGP breakup into the entropy poor HG occurs rapidly. 
The entropy conservation constraint fixes the final light 
quark yield. We assume a fast transition between QGP 
and HG phases, such that all hadron yields are at the 
same physical conditions as in QGP breakup. 
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Assuming that in the hadronization process the num- 
ber of b, c, s quark pairs doesn't change, the three un- 
known 7^, 7^, 7J 1 can be determine from their values in 
the QGP phase, 7Q, 7Q, 7^ (or dN? /dy) and the three 
flavor conservation equations, 



Table 1. Open charm, and bottom, hadron states we consid- 
ered. States in parenthesis either need confirmation or have 
not been observed experimentally, in which case we follow the 
values of Refs. [151130] . We implemented charm-bottom symme- 
try required for certain observables. Top section: D, B mesons, 
bottom section: D s , _B s -mesons 



dN? dN? 



dNj 
dy 



dy dy 
In order to conserve entropy: 

rfS H _ dS® 
dy dy 



s,c, 



dS_ 
dy ' 



(9) 



(10) 



a value 7^ ^ 1 is nearly always required when in the 

QGP phase j® G = 1 . This implies that yields of hadrons 
with light quark content are, in general, not in chemi- 
cal equilibrium, unless there is some extraordinary cir- 
cumstance allowing a prolonged period of time in which 
hadron reactions can occur after hadronization. Chemi- 
cal non-equilibrium thus will influence the yields of heavy 
flavored particles in final state as we shall discuss in this 
work. 

As noted at the beginning of this section, the use of the 
hadron phase space (denoted by H above) does not imply 
the presence of a real physical 'hadron matter' phase: the 
SHM particle yields will be attained solely on the basis of 
availability of this phase space as noted at the beginning of 
this section. Another way to argue this is to imagine a pot 
of quark matter with hadrons evaporating . Which kind 
of hadron emerges and at which momentum is entirely 
determined by the access to the phase space, and there 
are only free-streaming particles in the final state. 

Thinking in these terms, one can imagine that espe- 
cially for heavy quark hadrons some particles are pre- 
formed in the deconfined plasma, and thus the heavy hadron 
yields may be based on a value of temperature which is 
higher than the global value expected for other hadrons. 
For this reason we will study in this work a range 140 < 
T < 260 MeV and also consider sensitivity to this type 
of two-temperature chemical freeze-out of certain heavy 
hadron yield ratios. 



3 Relative charmed hadron yields 

3.1 Determination of js/jq 

We have seen considering s/S and also s and S individu- 
ally across the phase limit that in general one would ex- 
pect chemical non-equilibrium in hadronization of chemi- 
cally equilibrated QGP. We first show that this result mat- 
ters for the relative charm meson yield ratio D s / D, where 
D s (cs) comprises all mesons of type (cs) which are listed 
in the bottom section of table[IJ and D(cq) comprise yields 
of all (cq)states listed in the top section of table [TJ This 
ratio is formed based on the assumption that on the time 
scale of strong interactions the family of strange-charmed 
mesons can be distinguished from the family non-strange 
charmed mesons. 



hadron 


M[GeV] 


Q:c,b 


hadron 


M[GeV] 


g 


D u (0") 


1.8646 


Qu 


B u (0~) 


5.279 


1 


D+(Q-) 
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Qd 


B+(Q-) 
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1 


D*°(l-) 


2.0067 


Qu 


B*°(l-) 
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3 


D*+(l~) 


2.0100 


Qd 


B*+(l-) 


5.325 


3 


D°(0+) 


2.352 


Qu 


B°(0+) 


5.697 


1 


D+(0+) 


2.403 


Qd 


B+(0+) 


5.697 


1 


Df(l+) 


2.4222 


Qu 


Bf(l+) 


5.720 


3 


Dr + (1 + ) 


2.4222 


Qd 




5.720 


3 


D* 2 °(2+) 


2.4589 


Qu 


B?(2~) 


(5.730) 


5 


D*+(2+) 


2.4590 


Qd 


B*+(2+) 


(5.730) 


5 


Dt(0~) 


1.9868 


Qs 


B s u (0-) 


5.3696 


1 


D*+(l~) 


2.112 


Qs 


Bf(l-) 


5.416 


3 


D*J(0 + ) 

D:j(i + ) 


2.317 
2.4593 
2.573 


Qs 
Qs 
Qs 


B:°j(0+) 
Bfj(l + ) 
Bfj(2+) 


(5.716) 
(5.760) 
(5.850) 


1 

3 
5 



1.5 



T=0.14 GeV 
T=0.16 GeV 
T=0.18 GeV 



CO 






0.5 



Fig. 1. (Color on line) D/D a ratio as a function of 7^/7^ for 
T = 140 MeV (blue, solid line), T = 160 MeV (green, dashed 
line) and T = 180 MeV (red, dash-dot line). 



The yield ratio D s /D calculated using Eq. (jU and Eq. (O 
is shown in figure [1] Using Eq. (O we see that this ratio is 
proportional to 7JV7J 1 and weakly dependent on T: 



D 



7" £igpsimp m exp(~m Dsi /T) 
if S i g Di m i lf exp(~m Di /T) 



(11) 



A deviation of "f s /^/q from unity in the range we will see 
in section 15.41 leads to a noticeable difference in the ratio 
D s /D. We show in figure [T] results for T = 140, 160, 180 
MeV. In this T-domain, the effect due to lf/Y q mH ^ 1 is 
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Fig. 3. (Color on line) Jfi 4>(To) / D S D S ratio is evaluated at 
the dominant contribution to the variation of this relative two temperatures, T for heavy flavor hadrons, and Tb for (f> as 
yield. a function of T - To, with three values of T = 140, 160, 180 

MeV is considered with S H = S Q . 



3.2 Check of statistical hadronization model 

We next construct a heavy flavor particle ratio that de- 
pends on hadronization temperature only. To cancel the 
fugacities and the volume we consider the ratio J/&(f>/ D s D t 
in figure[2] Here J/\P yield includes the yield of (cc) mesons 
decaying into the J/& '. All phase space occupancies cancel 
since J/<P oc j^ 2 , <j) oc jf 2 , D s oc 7^7^ and similarly 
Ds oc 7^7^. When using here the particle D s {cs) and 
antiparticle D s (cs) any chemical potentials present are 
canceled as well. However, for the LHC and even RHIC 
environments this refinement is immaterial. 

This ratio J/^<j)/ D S D S , turns out to be practically con- 
stant, within a rather wide range of hadronization tem- 
perature T, see figure [2] The temperature range we study 
140 < T < 280 MeV allows us to consider an early freeze- 
out of different hadrons. To be sure of the temperature 
independence of J/&4>/ D S D S we next consider the possibil- 
ity that hadronization temperature T of charmed hadrons 
is higher than hadronization temperature Tq of <f>. We 
study this question by exploring the sensitivity of the ratio 
J/'Pcf)/ D S D S to the two temperature freeze-out in figure^ 
see bottom three lines for T = 180, 160, 140 MeV with j q 
from condition = S H , see figure HI If charmed hadrons 
hadronize later, T > Tq, and T — To < 60 the change in 
J/&(j)/ D S D S ratio is small (about 20%). If this were to be 
measured as experimental result, 



D S D„ 



(12) 



one could not but conclude that all particles involved are 
formed by mechanism of statistical hadronization. 



This interesting result can be understood, considering 
the behavior of the 7^(To)/7^(T) ratio, which increases 
rapidly with increasing T — Tq (see the top three lines in 
figure [3]). This ratio almost compensates the change in <fh 
yield, an effect we already encountered in the context of 
the results we show below in figure [5J For large T — Tq the 
ratio J/&(j)/D s D s begins to decrease more rapidly because 
j s increases for S H = , see figure ([6]). 



4 Entropy in Hadronization 
4.1 Entropy in QGP fireball 

The entropy content is seen in the final state multiplic- 
ity of particles produced after hadronization. More specif- 
ically, there is a relation between entropy and particles 
multiplicities, once we note that the entropy per particle 
in a gas is: 



^ = 3.61, 



^1=4, ^=4.2, (13) 
N ' N y ' 



for massless Bose, classical (Boltzmann) and Fermi gases, 
respectively. Effectively, for QGP with u, s, d, G degrees 
of freedom, S^/N® ~ 4 is applicable for large range of 
masses. Thus: 

-w* iJ w- (14) 

This in turn means that final state particle multiplicity 
provides us with information about the primary entropy 
content generated in the initial state of the QGP phase. 
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It is today generally believed that there is entropy con- 
serving hydrodynamic expansion of the QGP liquid. En- 
tropy is conserved in the fireball, and the conservation of 
entropy density a flow is expressed by: 



0, 



(15) 



where u^ 1 is local four velocity vector. A special case of 
interest is the so-called Bj0rken scenario [TB] for which 
Eq. ([15)) can be solved exactly assuming as initial con- 
dition scaling of the physical properties as a function of 
rapidity. This implies that there is no preferred frame of 
reference, a situation expected in very high energy colli- 
sions. Even if highly idealized, this simple reaction picture 
allows a good estimate of many physical features. Of rele- 
vance here is that the exact solution of hydrodynamics in 
(1+1) dimensions implies 



dS „ 
— = Const. 
dy 



(16) 



Thus entropy S is not only conserved globally in the 
hydrodynamic expansion, but also per unit of rapidity. 
Though we have (1+3) expansion, Eq. (fl6|) holds as long 
as there is, in rapidity, a flat plateau of particles yields. 
Namely, each of the domains of rapidity is equivalent, ex- 
cluding the projectile-target domains. However, at RHIC 
and LHC energies these are causally disconnected from 
the central rapidity bin, where we study the evolution of 
heavy flavor. The entropy we observe in the final hadron 
state has been to a large extent produced after the heavy 
flavor is produced, during the initial parton thermalization 
phase, but before strangeness has been produced. In or- 
der to model production of hadrons for different chemical 
freeze-out scenarios of the same reaction, we need to relate 
the entropy content, temperature and volume of the QGP 
domain. We consider for a u, d, G-chemically equilibrated 
QGP, and allowing for partial chemical equilibration of 
strangeness, the entropy content. 

The entropy density a can be obtained from the equa- 
tion 

ffS £ = -I^. (17) 



V 



V dT 



where the thermodynamic potential is: 

F Q (T, X q ,V) = -TlnZ( Q T, \ q ,V) ( 



(18) 



Inside QGP the partition function is a product of partition 
function of gluons Z g , light quarks Z q and strange quarks 
Z s , hence: 



Here g g is degeneracy factor for gluons and g q is degener- 
acy factor for quarks. The factor 7/4 = 2- 7/8 accounts for 
the difference in statistics and presence of both quarks and 
antiquarks. The number of degrees of freedom of quarks 
and gluons is influenced by strongly interactions, charac- 
terized by strong coupling constant a s : 



9g 



1 15 

1 - — -a s 

47T 



g q = 2 S 3 C 2/ (1- — a s 



50 



(22) 
(23) 



The case of strange quarks is somewhat more com- 
plicated, since we have to consider the mass, the degree 
of chemical equilibration, and guess-estimate the strength 
of QCD perturbative interactions. We have in Boltzman 
approximation: 



\nZ s = 2 p/a ^VT 3 , (24) 
g s = 2 s 3 c7s Q Q.5W(m s /T) (l - k— ) . (25) 

W(m/T) is function seen in Eq. ((6|). We allow both for 
strange and antistrange quarks, factor 2 p / a (which is for 
massless fermions 2-7/8 = 7/4). k at this point is a tem- 
perature dependent parameter. Even in the lowest order 
perturbation theory it has not been evaluated for massive 
quarks at finite temperature. We know that for massless 
quarks k ~ 2. Considering expansion in m/T, for large 
masses the correction reverses sign j!7| , which result sup- 
ports the reduction in value of k for m ~ T. We will use 
here the value k = 1 [T2"] . 

4.2 Number of degrees of freedom in QGP 

The entropy density following from Eq. (|17p is: 



4^ 
~90" 



7 



(9 a + - A 9 q )T A + —2 p /,g s T* 



A 
T ' 



(26) 



For strange quarks in the second term in Eq. (|26p we set 
the entropy per strange quarks to 4 units. In choosing 
S s /N s = 4 irrespective of the effect of interaction and 
mass value m s /T we are minimizing the influence of un- 
known QCD interaction effect. 

The last term in Eq. (|26|) comes from differentiation of 
the strong coupling constant a s in the partition function 
with respect to T, see Eq. ([17]) . Up to two loops in the (3- 
function of the renormalization group the correction term 
is US]: 



In Z = In Z g + In Z q + In Z s ; 
where for massless particles with X q = 1 



In Z„ 



90 
4 90 



VT 



(19) 

(20) 
(21) 



A = {b a a 2 s + bia 3 s ) 



2tt t 4 



rif5n r 



(27) 



with n f being the number of active fermions in the quark 
loop, rif ~ 2.5, and 



1 ( 2 



* = ^ (51 -£»/)• (28) 
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For the strong coupling constant a s we use 



a s {T c ) 



1 + C\xl{T/T c )' 



C = 0.760 ±0.002, 



(29) 



where T c = 0.16 GeV. This expression arises from renor- 
malization group running of a s {p), the energy scale at 
/i = 2ttT, and the value a s (Mz) = 0.118. A much more 
sophisticated study of the entropy in the QGP phase is 
possible [19] , what we use here is an effective model which 
agrees with the lattice data 20J . 

Eq. ([26|) suggests that we introduce an effective degen- 
eracy of the QGP based on the expression we use for en- 
tropy: 



9%{T)^g g {T) + -g q {T) + 2g s ^ 



A 90 



Which allows us to write: 



4^ 
~90~ 



*eff 



T 



and 



dS 
dy 



4tt2 
90 ; 



&T»^ c Const. 



(30) 



(31) 



(32) 



We show the QGP degeneracy in figure 3] as a func- 
tion of T e [140,260] MeV, top frame for fixed s/S = 
0, 0.03, 0.04 (from bottom to top), and in the bottom frame 
for the strangeness chemical equilibrium, J s = 1 (dashed) 
and approach to chemical equilibrium cases (solid) . When 
we fix the specific strangeness content s/S in the plasma 
comparing different temperatures we find that in all cases 
g® s increases with T. For s/S — we have a 2- flavor sys- 
tem (dotted line, red) and the effective number of degrees 
of freedom g® s varies between 22 and 26. The solid line 
with dots (green) is for s/S = 0.03, and the dot-dashed 
line (blue) gives the result for s/S = 0.04. 

In the bottom panel of figure [4] we note that like for 2 
flavors, case (s=0), for the 2+1-flavor system (j s = 1) g% 
increases with T (dashed line, black). g® s varies between 
30 and 35.5. The thin dashed lines indicate the range of 
uncertainty due to mass of the strange quark, which in 
this calculation is fixed with upper curve corresponding 
to m s — 90 MeV, and lower one m s = 160 MeV. The 
expected decrease in value of m s with T will thus have 
the effect to steepen the rise in the degrees of freedom 
with T . 

We now explore in a QGP phase the effect of an in- 
creasing strangeness fugacity with decreasing tempera- 
ture. This study is a bit different from the rest of this pa- 
per, where we consider for comparison purposes hadroniza- 
tion for a range of temperatures but at a fixed value of 
s/S. A variable jp(T) implies a more sophisticated, and 
thus more model dependent picture of plasma evolution. 
However, this offers us an important insight about g?L. 

We consider the function: 



V Q 



300 - T[MeV] 
160 ' 



(33) 




o a> 



32 



V 1 

m =125±35MeV 

s 

y s =(300-T)/160 
m =125±35 MeV 




160 



T [MeV] 



220 



Fig. 4. (color on-line) The Stefan-Boltzman degrees of free- 
dom g g based on entropy content of QGP, as function of 
temperature T. Upper frame: fixed s/S, the solid line with 
dots (green) is for a system with fixed strangeness per entropy 
s/S = 0.03, while dot-dashed (blue) line is for s/S = 0.04. The 
dotted (red) line is for 2-flavor QCD s/S = (u, d, G only); The 
bottom frame shows dashed (black) line 2+1-flavor QCD with 
m s — 125 ± 35 MeV (chemically equilibrated u,d,s,G system) . 
The (thick, thin) solid lines are for QGP in which strangeness 
contents is increasing as a function of temperature, see text. 



This is consistent with the kinetic computation of strange- 
ness production [IS]. At T = 140 MeV we have chemical 
equilibrium in the QGP phase, while and at the temper- 
ature T = 260 MeV we have 7 S = 0.25. The result for g% 
is shown as a thick (black) line in figure [31 with the range 
showing strange quark mass range m s — 125 ± 35 MeV. 
We see that in a wide range of temperatures we have 
29.5 < g% < 30.5. 

The lesson is that with the growth of 77 with decreas- 
ing T the entropy of the QGP is well described by a con- 
stant value gf s = 30 ± 0.5. Since the entropy is (nearly) 
conserved and g® s is (nearly) constant, Eq. (|32[) implies 
that we can scale the system properties using the con- 
straint T 3 dV/dy =Const. We stress again that these re- 
sults arises in a realistic QGP with 2 + 7^-flavors, but are 
model dependent and of course rely on the lattice mo- 
tivated description of the behavior of QGP properties. 
On the other hand it is not surprising that the rise of 
strangeness chemical saturation with decreasing tempera- 
ture compensates the 'freezing' of the q, G-degrees of free- 
dom with decreasing temperature. 
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Table 2. Reference values of volume, temperature, entropy, 
particle multiplicity 



dV/dy [im 3 ] 


T[MeV] 


dS Li /dy 


dN^/dy 


dN a /dy 


800 


200 


10,970 


2,700 


5,000 


2300 


140 


10,890 


2,700 


4,500 



4.3 Entropy content and chemical (non-)equilibrium 

We use as a reference a QGP state with dV/dy — 800 fm 3 
at T = 200 MeV, see table H We find from Eq. (0 the 
Q and H phase particle multiplicity. The hadron multi- 
plicity stated is what results after secondary resonance 
decays. The total hadron multiplicity after hadronization 
and resonance decays was calculated using on-line SHARE 
2.1 [21]. If a greater (smaller) yield of final state hadrons 
is observed at LHC, the value of dV/dy need to be revised 
up (down). In general expansion before hadronization will 
not alter dS/dy. We can expect that as T decreases, U 1 / 3 
increases. Stretching the validity of Eq. ([32jl to low tem- 
perature T = 140 MeV, we see the result in the second 
line of table [2] 

For QGP, in general the entropy content is higher than 
in a comparable volume of chemically equilibrated hadron 
matter, because of the liberation of color degrees of free- 
dom in the color-deconfined phase. The total entropy has 
to be conserved during transition between QGP and HG 
phases, and thus after hadronization, the excess of entropy 
is observed in excess particle multiplicity, which can be in- 
terpreted as a signature of deconfinement [22125] . The dy- 
namics of the transformation of QGP into HG determines 
how this additional entropy manifests itself. 

The comparison of entropy in both phases is temper- 
ature dependent but in the domain of interest i.e. 140 < 
T < 180 MeV the entropy density follows: 



3a* 



(34) 



Since the total entropy S is conserved or slightly increases, 
in the hadronization process some key parameter must 
grow in the hadronization process. There are two options: 
a) either the volume changes: 



3V 



h > 



V Q. 



(35) 



or 



b) the phase occupancies change, and since rii oc 7, 
q, s in hadron phase 



2.3 



7 H 



V3, jfhf>l. 



(36) 



In a slow, on hadronic time scale, transition, such as is 
the case in the early Universe, we can expect that case 
a) prevails. In high energy heavy ion collisions, there is 
no evidence in the experimental results for the long co- 
existence of hadron and quark phases which is required 
for volume growth. Consequently, we have U H ~ V® and 
a large value of 7J 1 is required to conserve entropy. The 
value of 7J 1 is restricted by 



7, cr =exp(m°/2T). 



(37) 



This value 7" is near to maximum allowed value, which 
arises at condition of Bose-Einstein condensation of pions. 
We will discuss quantitative results for 7^ (and 7J 1 ) below 
in subsection 15.41 



5 Strangeness in Hadronization 
5.1 Abundance in QGP and HG 

The efficiency of strangeness production depends on en- 
ergy and collision centrality of heavy ions collisions. The 
increase, with value of centrality (participant number), of 
per-baryon specific strangeness yield indicates presents of 
strangeness production mechanism acting beyond the first 
collision dynamics. The thermal gluon fusion to strangeness 
can explain this behavior [5] , and a model of the flow dy- 
namics at RHIC and LHC suggests that the QGP ap- 
proaches chemical equilibrium but also can exceed it at 
time of hadronization [12] . 

The strangeness yield in chemically equilibrated QGP 
is usually described as an ideal Boltzman gas. However, a 
significant correction is expected due to perturbative QCD 
effects. We implement this correction based on comments 
below Eq. (|25p. We use here the expression: 



dy 



-1? 



V it J dy 



(38) 



with the Boltzman limit density Eq. ([5]), and mass m s = 
125 MeV, g s = 6, A s = 1. The QCD correction corresponds 
to discussion of entropy in subsection 14.11 

We obtain strange quarks phase space occupancy 7J 1 
as a function of temperature from condition of equality 
of the number of strange quark and antiquark pairs in 
QGP and HG. Specifically, in the sudden QGP hadroniza- 
tion, quarks recombine and we expect that the strangeness 
content does not significantly change. For heavier flavors 
across the phase boundary this condition Eq. ([9]) is very 
well satisfied, for strangeness the fragmentation effect adds 
somewhat to the yield, 



dNl^dNl 
dy ~ dy 



(39) 



Using the equality of yields we underestimate slightly the 
value of strangeness occupancy that results. We recall that 
we also conserver entropy Eq. (| 10[) which like strangeness 
can in principle grow in hadronization, 



(40) 



using Eq. (JTUJ) we underestimate the value of 7^ 2 . 



Counting all strange particles, the number of pairs is: 



Ml 
dy 



H2 nt q ) 



^H 2(27 H^ q + n e q + Psn e q) 



(41) 
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0.035 



0.03 



CO 



0.025 



0.02 




140 



160 



180 



200 



220 



240 
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T[MeV] 



Fig. 5. (Color on line) Strangeness to entropy ratio s/S as 
function of temperature T, for the QGP (green, solid line for 
m s = 160 MeV, blue dash-dot line for m s = 90 MeV) with 
k = 1, see Eq. (|38[1 ; and for HG (light blue,dashed line) phases 
for 7q = 7 S = X q — A s = 1 in both phases. 



where n- q are densities of strange hadrons (mesons and 
baryons) calculated using Eq. |4]) in chemical equilibrium. 
P s is the strangeness content of the r\. The way we count 
hadrons is to follow strangeness content, for example — 



cq 



l K o 



cq 
l K- 



+ n c - q 



We impose in our calculations 
s — s. The pattern of this calculation follows an estab- 
lished approach, SHARE 2.1 [21] was used in detailed 
evaluation. 



5.2 Strangeness per entropy s/S 

Considering that both strangeness, and entropy, are con- 
served in the hadronization process, a convenient variable 
to consider as fixed in the hadronization process, is the 
ratio of these conserved quantities s/S. In chemical equi- 
librium we expect that in general such a ratio must be 
different for different phases of matter from which parti- 
cles are produced [24l[25lfT2] . 

We compare QGP and HG specific per entropy strange- 
ness content in figure O We show as function of temper- 
ature T the s/S ratios for chemically equilibrated QGP 
and HG phase. For the QGP the entropy S in QGP is 
calculated as described in section [4] and we use k = 1 in 
Eq. (|38|) . The shaded area shows the range of masses of 
strange quarks, considered, results for m s — 90 MeV (up- 
per (blue) dash-dotted line) and m s — 160 MeV ((green) 
solid line) form the boundaries. The central QGP value is 
at about s/S = 0.032. 

The short-dashed (light blue) line shows the hadron 
phase s/S value found using SHARE 2.1 . For HG near to 
usual range of hadronization temperature T ~ 160 MeV 
we find s/S ~ 0.025. In general formation of QGP implies 



and increase by 30% in s/S. Both HG and QGP phases 
have a similar specific strangeness content at T = 240- 
260 MeV, however it is not believed that a HG at such 
high temperature would be a stable form of matter. This 
HG to QGP dissociation, or QGP hadronization depends 
on the degree of strangeness equilibration in plasma [29J , 
and other dynamical factors. 

In the QGP the value of s/S for the range of realistic 
hadronization temperature 140 < T < 180 MeV is in gen- 
eral larger than in HG. This implies that generally, the 
abundance of strange hadrons produced in hadronization 
over saturates the strange hadron phase space, if QGP 
state had reached (near) chemical equilibrium. Moreover, 
since we are considering the ratio s/S and find in QGP 
a value greater than in HG, for chemical equilibrium in 
QGP the hadronization process will lead to 7^/7^ > 1. 

One can wonder if we have not overlooked some dy- 
namical or microscopic effect which could adjust the value 
of s/S implied by QGP to the value expected in HG. 
First we note that the fast growth of the volume V cannot 
change s/S. Moreover, any additional strangeness produc- 
tion in hadronization would enhance the over-abundance 
recorded in the resulting HG. Only a highly significant en- 
tropy production at fixed strangeness yield in the hadroniza- 
tion process could bring the QGP s/S ratio down, mask- 
ing strangeness over-saturation. A mechanism for such en- 
tropy production in hadronization is unknown, and more- 
over, this would further entail an unexpected and high 
hadron multiplicity excess. 

One could of course argue that the perturbative QCD 
properties in the QGP are meaningless, the entropy in 
QGP is much higher at given temperature. However, the 
properties of QGP have been checked against the lattice 
results, and the use of lowest order expressions is justified 
in these terms [20]. Moreover, the value of s/S is estab- 
lished way before hadronization. 



5.3 Wroblewski ratio W s 

At this point it is appropriate to look at another observ- 
able proposed to study strangeness yield, the Wroblewski 
ratio [26] : 

(42) 



+ (dd) ' 



W s compares the number of newly produced strange quarks 
to the produced number of light quarks. In an equili- 
brated deconfined phase W s compares the number of ac- 
tive strange quark degrees of freedom to the number of 
light quark degrees of freedom. 

The ratio s/S compares the strange quark degrees of 
freedom to all degrees of freedom available in QGP. There- 
fore as function of T the ratios s/S and W s can behave 
differently: Considering the limit T — > T c a constant s/S 
indicates that the reduction of s-degrees of freedom goes 
hand in hand with the 'freezing' of gluon degrees of free- 
dom, which precedes the 'freezing' of light quarks. This 
also implies that for T — > T c in general W s diminishes. 



10 Inga Kuznetsova and Johann Rafelski: Heavy Flavor Hadrons in Statistical Hadronization of Strangeness-rich QGP 



The magnitude of m s , the strange quark mass decisively 
enters the limit T — * T c . 

For T » T c the ratio W s can be evaluated comparing 
the rates of production of light and strange quarks, using 
the fluctuation-dissipation theorem [27 , which allows to 
relate rate of quark production to quark susceptibilities 
Xi (see Eqs. (11) and (12) in [27]): 



W s ~IL 



2Xs 



Xu+Xd 



(43) 



An evaluation of R x as function of temperature in lat- 
tice QCD has been achieved [25]. For T ~ 2.5T C the result 
obtained, W s — > R x ~ 0.8, is in agreement with the expec- 
tation for equilibrium QGP with nearly free quarks, with 
mass of strangeness having a small but noticeable signifi- 
cance. With decreasing T, the ratio R x — > 0.3 for m s = T c . 
However, this value of m s is too large, the physical value 
should be nearly half as large, which would result in a 
greater R x . Moreover, for T — > T c the relationship of W s 
to R x , Eq. (f43|) is in question in that the greatly reduced 
rate of production of strangeness may not be satisfying 
the conditions required in Ref. [27] . 

Comparing the observables s/S and W s we note that 
the experimental measurement requires in both cases a de- 
tailed analysis of all particles produced. At lower reaction 
energies there is additional complication in evaluation of 
W s due to the need to subtract the effect of quarks brought 
into the reaction region. Turning to the theoretical com- 
putation of s/S and W s we note that the thermal lattice 
QCD evaluation of s/S is possible without any approxi- 
mation, even if the actual computation of entropy near the 
phase boundary is a challenging task. On the other hand, 
the lattice computation of W s relies on production rate 
of strangeness being sufficiently fast, which cannot be ex- 
pected near to the phase boundary. Moreover, the variable 
s/S probes all QGP degrees of freedom, while W s probes 
only quark degrees of freedom. We thus conclude that s/S 
is both more accessible theoretically and experimentally, 
and perhaps more QGP related observable, as compared 
to W s , since it comprises the gluon degrees of freedom. 



5.4 Strangeness chemical non-equilibrium 

In order that in fast hadronization there is continuity 
of strangeness Eq. (f3"5]) . and entropy, Eq. ([4TJ]) the hadron 
phase -ff^l and 7^ 7^ 1. We have to solve for 7^ and 
7J 1 simultaneously Eqs. (|39I40|) . 

In figure [5] we show as a function of T the strange 
phase space occupancy 7^, obtained for several values of 
s/S ratio (from top to bottom 0.045, 0.04, 0.035, 0.03, 
0.025) evaluated for S" Q = S" H . The solid line shows 7" 
for = s H and S^ = S^. The maximum allowed value 
Eq. ([3"T]) is shown dashed (red) . 

In figure [7] we show results for 7^ /7J 1 (where 7J 1 = 1 
we show 7I 1 ) . We consider the three cases: j q — 1 , 7J 1 = 
7^ r , and entropy conservation S H — S® for s/S — 0.045, 
s/S = 0.04, s/S = 0.035, s/S = 0.03, s/S = 0.025 (dash- 
dot lines) (lines from top to bottom). We see that except in 
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6. (color on-line) Phase space occupancy as a function 

H AO„„ „„iiJ 1:„„\ „,H 



Fig 

of T: (blue, solid line), 7^ (dash-dotted lines, from top to 
bottom) for s/S = 0.045, for s/S = 0.04 (thick line), s/S = 
0.035, s/S = 0.03 (thick line), s/S = 0.025; 7" (red, dashed 
line) . 



case that strangeness were to remain well below chemical 
equilibrium in QGP (s/S ~ 0.03), the abundance of heavy 
flavor hadrons we turn to momentarily will be marked by 
an overabundance of strangeness, since practically in all 
realistic conditions we find 7^ > 7^. 

In figured] we show s/S ratio as function of 7^/7^- 
The solid line is for T = 200 MeV, 7" = 0.83, S Q = S H , 
dashed line for T = 170 MeV, 7" = 1.15 S Q = S H and 
dash-dot line for T = 140, 7J 1 = 1.6 MeV, S Q = S H . We 
also consider j q — 1 case for T — 170 MeV (dot marked 
(purple) solid line). In this case strangeness content 7 s /7 g 
is higher than for S H = S® with the same T and s/S. 
In the limit = 7^, Eq. [37] (T = 200 MeV, solid, thin 
line; T = 170 MeV, dashed line) the strangeness content 
7?/7g i i s minimal for given T and s/S. 

These results suggest that it is possible to measure the 
value of s/S irrespective of what the hadronization tem- 
perature may be, as long as the main yield dependence is 
on the ratio 7^/7^- Indeed, we find that the ratio (f>/K + ; 



Is 



eq 



K+ 



vH cq 



(44) 



is less sensitive to hadronization temperature compared to 
its strong dependence on the value of s/S. In figure [5] we 
show the total hadron phase space ratio <f>/K + as function 
of T for several s/S ratios, and for 7^ = 1 (chemical 
equilibrium, dashed (red) line). The K + yield contains 
the contribution from the decay of <fi into kaons which is 
a noticeable correction. 
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Fig. 7. (color on line) 7^/7? (=7? at 7J 1 = 1) as a function 
of hadronization temperature T. Top frame: = 1, middle 
frame: = 7"", and bottom frame: S H = S®. Lines, from top 
to bottom: s/S = 0.04 (blue, solid line), s/S — 0.03 (green, 
dashed line), s/S — 0.025 (red, dash-dot line) 



Table 3. Specific and absolute strangeness yield for different 
reaction volumes at T = 200. 



s/S 


ds/dy 


dV/dy [im-' 6 } 


T [MeV] 


0.045 


550 


1000 


200 


0.04 


360 


800 


200 


0.035 


250 


700 


200 


0.03 


165 


600 


200 


0.025 


106 


500 


200 


0.022 


83 


500 


200 



We record in table [3] for given s/S and volume dV/dy 
the corresponding total yields of strangeness, which may 
be a useful guide in consideration of the consistency of 
experimental results with what we find exploring heavy 
flavor hadron abundance. 



0.06 1 1 1 1 1 1 1 1 1 1 1 1 r 




Fig. 8. (color on line) Strangeness to entropy ratio, s/S, as 
a function of j s /j g for T = 200 MeV, 7, = 0.083, S H = 
S Q (solid line), T = 170 MeV, j q = 1.15, S H = S Q (dashed 
line), T = 140 MeV, 7, = 1.6, S H = S Q (dash-dotted line); 
7 9 = 1 (dot marked solid); -y q = 7, r : T = 200 MeV (thin solid 
line), T = 170 MeV (thin dashed line). 
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Fig. 9. (Color on line) The ratio 4>/K + as a function of T. 
Dashed line (red) is for chemical equilibrium. Solid line with 
dots (green) s/S = 0.03, solid line (blue) s/S = 0.04, dash-dot 
line (per) is for s/S = 0.022. 
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6 Yields of heavy flavored hadrons 

6.1 Phase space occupancy 7^ and 7^ 



The first step in order to determine the yields of heavy 
flavor hadronic particles is the determination of the phase 
space occupancy 7^ and 7^. 7 
of number of these quarks (i.e 
pairs) in QGP and HG. The yield constraint is 



H is obtained from equality 
of quark and anti quark 



dN c dV r 11 



11 2 



(n ceq 



2 7 % 



eq 

ccq 



cq \ 

CCS J 



where open 'op' charm yield is: 



(45) 



H eg 



n op= lq n D+ls 



n Ds + 7 



H2 , 
9 



'qqc 



,H/ n oq 
lq sqc 



Here and 



cq 



are densities of D and D< 



"b 7s n ssc- 

(46) 
mesons, 



respectively, in chemical equilibrium, n^ c is equilibrium 
density of baryons with one charm and two light quarks, 
nf^ c is density of baryons with one charm (or later on 
one bottom quark) and two strange quarks (J?°, ■!?£) in 
chemical equilibrium and n^? d is equilibrium particle den- 
sity with both, a charm (or bottom) and an anticharm 
(or antibottom) quark (C=0, B=0, S=0). The equilibrium 
densities can be calculated using Eq.Q. 7* 1 can now be 
obtained from Eq. ([45|) . 

Similar calculations can be done for 7J 1 . The only dif- 
ference is that we need to add number of B c mesons to 
the right hand side of Eq. (|45|) , 



(47) 



dN Bc h h cq dV 



cq 
*J3c 



is density in chemical equilibrium of B c . In the cal- 
culation of 7<P the contribution of term with ub c is very 
small and we did not consider it above. 

The value of j c is in essence controlled by the open sin- 
gle charm mesons and baryons. For this reason we do not 
consider the effect of exact charm conservation. The rela- 
tively small effects due to canonical phase space of charm 
are leading to a slight up-renormalization of the value of 
7 C so that the primary dN c /dy yield is preserved. This 
effect enters into the yields of multi-charmed and hidden 
charm hadrons, where the compensation is not exact and 
there remains slight change in these yields. However, the 
error made considering the high yield of charm is not im- 
portant. On the other hand for multi-bottom and hidden 
bottom hadrons the canonic effect can be large, depend- 
ing on actual bottom yield, and thus we will not discuss 
in this paper yields of these hadrons, pending extension 
of the methods here developed to include canonical phase 
space effect. 

We consider in figure [TU] the temperature dependence 
of both 7^ (top) and 7J 1 (bottom) for the heavy flavor 
yield given in Eqs. (|1I2[) . In the non-equilibrium case (solid 
lines) the space occupancy 7^ is obtained from Eq. (I41|) 
and 7^ is chosen to keep S H — S®. 7^ depend on 
7 S and j q : the value of N s in Eq. (j4"Tj) is chosen to have 
s/S = 0.04 after hadronization, the corresponding 7^ and 
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Fig. 10. (Color on line) 7{T(b = 1) (upper panel), and 7 J 1 
(c=10) (lower panel) , as functions of temperature of hadroniza- 
tion T. The solid lines are non-equilibrium for s/S — 0.04 
with S® = S H , dashed lines are equilibrium case 7 S =7 9 = 1 
and dot-dash lines are for s/S = 0.04 with maximal value of 
7,(7, =77) ( d V/dy = 800 fm 3 ). 



7^ are shown in figures [5] and [7J. Since applicable 7^ may 
depend on hadronization dynamics and/or details of equa- 
tion of state of QGP, we show charm quark phase space oc- 
cupancies also for maximum possible value of 7J 1 — > 7^ r ), 
also considered at s/S = 0.04 for all hadronization tem- 
peratures. We can compare our results with the chemi- 
cal equilibrium (dashed lines) setting 7^ = 7J 1 = 1 in 
Eq. At hadronization condition T = 160 ± 20 MeV 
temperatures we see in figure [TOl a significant (considering 
the fast changing logarithmic scale) difference between the 
chemical equilibrium, and non-equilibrium (s/S=0.04) re- 
sults. 

In figure[TT]we show the ratio 1 UM) eq as a function 
of hadronization temperature T. This helps us understand 
when the presence of chemical nonequilibrium is most no- 
ticeable. This is especially the case should heavy flavor 
hadronization occur at the same temperature T = 140- 
170 MeV as is obtained for non-heavy hadrons, and/or 
when the entropy content of light hadrons is maximized 
with 7J 1 — > 7". When no additional entropy is formed 
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Fig. 11. (Color on line) llbeq an d 7c llceq-t as functions of 
temperature of heavy flavor hadronization T. The solid with 
dot marks line is for y^/y^g with s/S = 0.04, solid line is 
for 7,? with s/S = 0.04, dot-dash and dashed lines are for 
s/S = 0.04 with maximal value of y q — ► y q r for 7J 1 /y^ eq and 
for yf/yf eq , respectively. 



in hadronization, that is S H = S®, lfn,)/lf(h) eq exceeds 
unity for T > 200 MeV, at which point the heavy fla- 
vor hadron yields exceed the chemical equilibrium expec- 
tations. In general we find that heavy hadron yields if 
produced at normal hadronization temperature would be 
effectively suppressed, compared to statistical equilibrium 
results, by the high strangeness yield. This happens since 
the phase space is bigger at 7^ > 1, and thus a smaller 
y^ b is required to reach a given heavy flavor yield. 

7J 1 and 7J 1 are nearly proportional to dN^ c /dy, re- 
spectively. The deviation from the proportionality is due 
to the abundance of multi-heavy hadrons and it is small. 
To estimate this effect more quantitatively we first evalu- 
ate: 

dl% (dV 



7 c0 



(48) 



i.e. the value expected in absence of multi heavy hadrons. 
Next we compare with the result when we take into ac- 
count the last three terms in Eq. (|45[) . The influence of 
these therms depend not only on dN c /dy but also on 
dN c /dy/dV/dy. For fixed dV/dy = 800 fm -3 in the range 
of dN c /dy = (5,30), we find that y^/N c (and therefore 
yields of open charm hadrons) changes at temperature 
T = 140 MeV by - 6% for the s/S = 0.04. For the chem- 
ical equilibrium case 7 S = y q — 1, yjr/N c changes up to 
15% at the same conditions. For the particles with hidden 
charm or 2 charm quarks the yields are proportional 7?, 
therefore changes in their yields will be about twice larger. 
For RHIC N c < 3 and dV/dy = 600 fm" 3 the dependence 
of yields on N c is much smaller. 



The multiplicity dN c /dy can also influence 7^, since as 
we noted it also includes a term proportional to Jc n< B C - ^ n 
the range of N c — (5,30), yb/Nb changes at temperature 
T = 0.14 MeV by - 0.5% for s/S = 0.04. Since the mass 
of &-quark is much larger than that of c-quark, the effect 
due to multi-bottom states is negligible. 



6.2 D, Ds, B, Bs meson yields 

In next sections we will mostly consider particles yields 
after hadronization and we will omit superscript H in 
7s. Considering Eq. ([4]), we first obtain y c as a function 
of y s /y q ratio and T. Substituting this y c and appropri- 
ate equilibrium hadron densities into Eq. ^ we obtain 
yields of D(B) and D S (B S ), as functions of the y s /yq ra- 
tio, at fixed temperature, which are shown on figure [T2l In 
the upper panel we show the fractional yields of charmed 
D/N c and D 8 /N c mesons, and in the lower panel B/Nf, 
and B s /N b for T = 200 MeV (solid line), T = 170 MeV 
(dashed line), T = 140 MeV (dash-dot line). Fractional 
yield means that these yields are normalized by the to- 
tal number of charm quarks N c and, respectively bottom 
quarks Nb, and thus tell us how big a fraction of avail- 
able heavy flavor quarks binds to non-strange and strange 
heavy mesons, respectively. Using figure [5] the ratio 7 s /7 g 
can be related to the s/S ratio. y q was chosen to con- 
serve entropy during hadronization process, see figure [6] 
In general the heavy non-strange mesons yield decreases 
and strange heavy meson yield increases with 7 s /7 9 . The 
yields D, B and D s , B s are sum over exited states of D, B 
and D s , B s respectively, see table[T]for the 'vertical tower' 
of resonances we have included. 

Using 7c, 7 S , 7 g at a given T (see figures [H [71 [lOj) we 
have now all the inputs required to compute absolute and 
relative particle yields of all heavy hadrons which we can 
consider within the grand canonical phase space. When 
we consider chemical equilibrium case, we use naturally 

ls=lq = 1. 

In figure [13] we consider the yields shown in figure [12] 
as a functions of hadronization temperature. The dashed 
blue and green lines were obtained for chemical equilib- 
rium yields of D and D s respectively. The extreme upper 
and lower lines are for fractional D and D s yields with 
s/S — 0.03 (dot marked, blue and green lines, respec- 
tively), while the central lines are for s/S = 0.04 (solid, 
blue and green lines). Also we show fractional yields for 
maximal possible value y q — > r ) q r , see figure [6] for 7^ r (T)) 
(dash-dot lines) and Eq. (f37|) . 

We note that there is considerable symmetry at fixed 
T between the fractional yields of charmed, and bottom 
mesons, for the same condition of s/S. The chemical equi- 
librium results show significant difference between strange 
and non-strange heavy mesons. In the case of chemical 
equilibrium, for the considered very wide range of hadroniza- 
tion temperatures D s /N c ~ B s /Nb ~ 0.2 are nearly con- 
stant. A significant deviation from this result would sug- 
gest the presence of chemical non-equilibrium mechanisms 
of heavy flavor meson production. 
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Fig. 12. (Color on line) Upper panel, fractional charm meson 
yield, and lower panel, fractional bottom meson yields as a 
function of 7s /7 ? ratio for fixed hadronization temperature T. 
Upper lines in each panel are for D(B) mesons, solid line is 
for T = 200 MeV, -y q = 1.1, dashed line is for T = 170 MeV, 
7, = 1.15 and dash-dot line is for T = 140 MeV, j q = 0.83 
(S H = S Q ). 



The yields of D s /N c and B s /Nb are very similar, and 
similarly so for D/N c and B/Nj,. Thus the relative yield of 
either of these mesons measures the relative yield of charm 
to bottom participating in the statistical hadronization 
process: 

D, D N r 

(49) 



B„ 



B N b 



This is a very precise result, which somewhat depends 
on the tower of resonances included, and thus in particu- 
lar on the symmetry in the heavy quark spectra between 
charmed and bottom states which we imposed. 

It is useful to reconsider here the ratio D/D s (B/B s ) 
which is proportional to 7 g /7 s , see figure [T]for D s /D pre- 
sented as a function of ^s/lq- We consider this ratio now 
as a function of T, the upper panel in figure [14] is for 
charm, the lower for bottom. We see that there is consid- 
erable symmetry in the relative yields between charmed 
and bottom mesons with upper and lower panels looking 
quasi-identical. Except for accidental values of T where 
the equilibrium results (blue, dashed lines) cross the fixed 
s/S results, there is considerable deviation in these ra- 
tios expected from chemical equilibrium. For LHC with 
s/S = 0.04 this ratio is always noticeably smaller than 
in chemical equilibrium (solid purple line is for = S H 
and purple, dash-dot line is for 7 g = 7^ r ). Even for RHIC- 
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Fig. 13. (Color on line) Upper panel, fractional charm me- 
son yield, and lower panel, fractional bottom meson yields. 
Equilibrium (dashed lines) and non-equilibrium for s/S = 0.03 
(point marked solid line) and s/S = 0.04 (solid line) for D/N c 
(blue lines, upper panel); D s /N c (green lines, upper panel); for 
D/N c and D s /N c with s/S = 0.03 and s/S = 0.04 for 7 , = 7 = r 
(dash-dotted lines); B/Nb (solid line, lower panel); and B s /Nb 
(point marked solid line, lower panel), for B/Nt and B s /Ni, 
with s/S = 0.03 and s/S = 0.04 for -y q = -y c q r (dash-dot lines); 
as a function of T. 



like conditions with s/S = 0.03 this ratio is smaller than 
in chemical equilibrium for all temperatures when entropy 
conservation in hadronization is assumed, S® = S (dot 
marked solid, green line). 



6.3 Heavy baryon yields 

As was the case comparing charm to bottom mesons we 
also establish a symmetric set of charmed and bottom 
baryons, shown in the table[U Many of the bottom baryons 
are result of theoretical studies and we include that many 
states to be sure that both charm and bottom are con- 
sider in perfect symmetry to each other. In figure [15] (up- 
per panel) we show hadronization temperature dependen- 
cies of yields of baryons with one charm quark normalized 
to charm multiplicity N c . We show separately yields of 
baryons without strange quark (A c + S c )/N c , and with 
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Fig. 14. (Color on line) Ratios D/Ds (upper panel) and B/Bs 
(lower panel) are shown as a function of T for different s/S 
ratios and in chemical equilibrium. Solid line is for s/S — 0.04, 
dash-dot line is for s/S — 0.03, dashed line is for jf = jf = 1. 

Table 4. Charm and bottom baryon states considered. States 
in parenthesis are not known experimentally and have been 
adopted from theoretical source [31] . 



hadron 


M[GeV] 


Q : c,b 


hadron 


M[GeV] 


g 


/l c +(l/2+) 


2.285 


udQ 


A0(l/2+) 


5.624 


2 


At(l/2~) 


2.593 


udQ 


A0(l/2") 


(6.00) 


2 


A+(3/2-) 


2.627 


udQ 


A0(l/2") 


(6.00) 


2 


E+(l/2+) 


2.452 


qqQ 


Z° b (l/2+) 


(5.77) 


6 




2.519 


qqQ 


S° b *(3/2+) 


(5.78) 


12 


H c (l/2+) 


2.470 


qsQ 


3,(1/2+) 


(5.76) 


4 


s' c (l/2+) 


2.574 


qsQ 


s£ (1/2+) 


(5.90) 


4 


3 c (3/2+) 


2.645 


qsQ 


3^(3/2+) 


(5.90) 


8 


Qc{l/2 + ) 


2.700 


ssQ 


n„(i/2+) 


(6.00) 


2 


n c (3/2+) 


(2.70) 


ssQ 


0,(3/2+) 


(6.00) 


4 



one strange quark S=l (S c /N c ). We show two cases for 
s/S = 0.04 with conserved entropy at hadronization S® = 
5 H (solid lines) and with maximum possible entropy value 
"f q — (dash-dot lines). The chemical equilibrium case 
7 g = 7 S = 1 is also shown (dashed lines). The upper lines 
of each type are for (A c + S c )/N c , the lower lines are for 
S c /N c . A similar result is presented for bottom baryons 
in the lower panel of figure [15] We note that the result 
for bottom baryons is more uncertain since most baryon 
masses entering are not experimentally verified. 
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Fig. 

0.04, S Q = S H (solid lines), s/S = 0.04, 7, = 7= r , (the upper 
panel) upper lines for each type are for ratio (A c + E c )/N c 
and lower lines are for S c /N c (upper panel) and (lower panel) 
upper lines of each type are for (A b + S b )/N c and lower lines 
are for S b /N b as functions of T. 



We note that the results shown figure [15] imply that 
under LHC conditions at least 15% of heavy flavor can 
be bound in heavy baryons, but possibly 30%. For large 
7 g = jg T > 1 we see increase in (A c + E c )/N c yields com- 
pared to chemical equilibrium and especially compared 
to entropy conserved hadronization = S H . This is so 
since yields are proportional to 7*, 7s7 9 . This results to 
relative suppression the D s /N c (see figure [T3|). 

In figure [TBI we show ratio cqq/cqs = (A c + E c )/E c 
as a function of 7^/7, for T = 200 MeV (dash-dot line), 
T = 170 MeV (solid line) and T = 140 MeV (dashed 
line). This dependence is linear, the slope depends only 
on hadronization temperature T. The 7s /7g ratio can be 
converted to s/S ratio using figured] 

The yield of multi-strange charmed baryon, Q c (css) is, 
similar to the light multi-strange hadrons, much more sen- 
sitive to chemical non-equilibrium. In figure [T7] we see a 
large increase in fractional yield of fi c {css)/N c for s/S = 
0.04 and 5^ = 5 H (solid line) compared to the chemical 
equilibrium (dashed line) expectation for the entire con- 
sidered range of hadronization temperature. As expected, 
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Fig. 16. (Color on line) The ratios cqs/cqq = S C /(A C + E c ) 
(upper lines) and css/cqq = (2 C /(A C + S c ) (lower lines) for 
T = 200MeV (dash-dot line), T = 170 MeV (solid line) and 
T = 140 MeV (dashed line) as functions of j s /"fq ■ 



Table 5. Hidden charm and multi heavy hadron states con- 
sidered. States in parenthesis are not known experimentally 



hadron 




mass(GeV) 


g 


T) C (1S) 


CC 


2.9779 


1 


J/P(1S) 


CC 


3.0970 


3 


Xeo(lP) 


CC 


3.4152 


1 




CC 


3.5106 


3 


h c (lP) 


CC 


3.526 


3 


Xta(lP) 


CC 


3.5563 


5 


Vc(2S) 


CC 


3.638 


1 




CC 


3.686 


3 


i> 


CC 


3.770 


3 


Xc2(2P) 


CC 


3.929 


5 


i> 


CC 


4.040 


3 


i> 


CC 


4.159 


3 


i> 


CC 


4.415 


3 


B c 


be 


6.27 


1 


' — <cc 


ccq 


3.527 


4 


^cc 


ccs 


(3.660) 


2 



6.4 Yields of hadrons with two heavy quarks 
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Fig. 17. (Color on line) Q c (css)/N c as function of T: dashed 
line for chemical equilibrium; solid lines are for S Q = S s , 
dashed dotted lines are for 7 9 = 7"": both for s/S — 0.03 
and s/S = 0.04 (upper lines). 



this yields increase with T. This also means that higher 
formation temperature can be invoked to explain an un- 
usually high yield. We expect that at LHC more than one 
percent of total charm yield will be found in the £2 c (css) 
state. 



We consider multi- heavy hadrons listed in the table [5] The 
yields we will compute are now more model dependent 
since we cannot completely reduce the result, it either re- 
mains dependent on the reaction volume dV/dy, or on the 
total charm(bottom) yields dN/dy. For example the yields 
of hadrons with two heavy quarks are approximately pro- 
portional to 1/ (dV I dy) because for heavy quarks is 
proportional to 1/dV/dy, see Eq. (|45|) : 



(50) 
(51) 



dy ^ c dy dV/dy 1 
dN Bc H H dV 1 
dy K7 ° lb dy X dV/dy 



Moreover, unlike it is the case for single heavy hadrons, the 
canonical correction to grand-canonical phase space does 
not cancel out in these states, adding to the uncertainty. 

Thus the result we present must seen as a guiding the 
eye and demonstrating a principle. In figure [18] we show 
the yield of hidden charm cc mesons (see table [5]) normal- 
ized by the square of charm multiplicity N% as a function 
of hadronization temperature T. We consider again cases 
with s/S = 0.03 (upper panel) and s/S — 0.04 (lower 
panel), solid line is for S H — S®, dot-dash line is for 
7 g = 7^ r , and dot-dash line is for -f q — 7?". The chemical 
equilibrium cc mesons yields are shown (dashed lines on 
both panels) for two different values of dV/dy — 600 fm 3 
for T = 200 MeV (upper panel) and dV/dy = 800 fm 3 for 
T = 200 MeV (lower panel). 

The yield of cc mesons is much smaller for s/S = 0.04 
than in equilibrium for the same dV/dy for large range 
of hadronization temperatures. For s/S — 0.03 the ef- 
fect is similar, but suppression is not as pronounced. For 
7 g = 7g r suppression the yield of hidden charm particles 
is always smaller than equilibrium value. This suppression 
occurs due to competition with the yield of strange-heavy 
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Fig. 18. (Color on line) cc/N^ yields as a function of 
hadronization temperature T, at dV/dy — 600 fm -3 for T = 
200 MeV, s/S = 0.03 (upper panel), dV/dy = 800 fm" 3 for 
T = 200 MeV, s/S = 0.04 (lower panel). Results shown are for 
S® — S H (solid lines), for = j q r (dash-dot lines), and for 
chemical equilibrium case (dashed lines, s/S is not fixed). 



mesons, and also, when j q > 1, with heavy baryons with 
two light quarks. The enhanced yield of D, D s and heavy 
baryons in effect depletes the pool of available charmed 
quark pairs, and fewer hidden charm cc mesons are formed. 
For particles with two heavy quarks the effect is larger 
than for hadrons with one heavy quark and light quark(s). 

In figure [19] we compare the J ft yield to the chemical 
equilibrium yield \P/Jft eq , as a function of 7^/7^, each 
line is at a fixed value 7^. This ratio is: 



N-, 



hid 



Jft 

Jftcq Nhid eq 



/ c cq 



(52) 



Jft/Jft cq always decreases when js/lq increases. For j q = 
7 cr J ft I Jft C q is smaller than unity even when 7 S — * 0, 
because of large phase space occupancy of light quarks. 
Jft/Jft cq > 1 for small 7, and small 7 s /7 9 . 

Considering the product of J ft and <j> yields normal- 
ized by N 2 we eliminate nearly all the uncertainty about 
the yield of charm and/or hadronization volume. How- 
ever, we tacitly assume that both J ft and <f> hadronize at 
the same temperature. In figure |2"01 we show J/tyfi/N 2 as 
function of "f s /"f q . There is considerable difference to the 
ratio considered in figure [21 We see mainly dependence on 
7 s /7 g . As before, see section [X2l J^ is the sum of all states 
cc from table[5]that can decay to Jft. We show results for 
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Fig. 19. (Color on line) Ratio Jft/Jft eq — Jc/lceq as a func- 
tion of 7JV7J 1 at fixed value of 7^ and if required, entropy 
conservation. Shown are: T = 200 MeV at j q — 0.83 (dot- 
dash line) and at j q = "/ q r = 1.42 (lower solid line (purple) ); 
T — 170 MeV at -y q = l(upper dashed line) , at 7, = 1.15, 
(upper solid line (red)), and at 7, = 7, r = 1.51, (lower dashed 
line); and T = 140 MeV, j q = 1.6 



T = 200 MeV (solid lines), T = 170 MeV (dashed line) 
and T = 140 MeV (dash-dot line). The j q , for each T, is 
fixed by entropy conservation condition during hadroniza- 
tion (figure [6|) (thick lines) or by j q = 7^ r (thin lines) . 
For T = 140 MeV these lines coincide. T = 170 MeV, 
jq = 1 case is also shown (solid line with dot markers). 
The s/S values, which correspond to given js/jq ratio can 
be found in figured] Figure l20l shows that despite the yield 
4>/{dV/dy) increasing as {j s /lq) 2 , Jft(j)/N 2 is increasing 
as js/lq considering compensation effects. 

A similar situation, as in figure [19] for hidden charm, 
arises for the B c meson yield, see figure[2T] where B c /N c Nb 
ratio is shown as a function of hadronization tempera- 
ture T, for the same strangeness yield cases as discussed 
for the hidden charm meson yield. Despite suppression in 
strangeness rich environment, the B c meson yield contin- 
ues to be larger than the yield of B c produced in single NN 
collisions, where the scale yield is at the level of ~ 10 -5 , 
see cross sections for bb and B c production in [llj and in 
[32j . respectively. 

In figure [H] we show N 2 scaled yields of ccq and ccs 
baryons as a function of temperature. Upper panel shows 
aside of the equilibrium case (dashed lines) the yields for 
s/S = 0.04 with S H = S Q (solid lines) and with lq = 7 = r 
(dash-dot line) for dV/dy = 800 fm~ 3 for T = 200 MeV. 
Lower panel is for dV/dy = 600 fm -3 and s/S = 0.03. For 
the ccg baryons the chemical non-equilibrium suppression 
effect is similar to what we saw for cc and Bc mesons. 
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Fig. 20. (Color on line) J/9(j)/N^ states yields as a function 
of 7 S /7 9 ratio for T = 200 MeV, S Q = S H (solid line) and 
7„ = 7j r (solid thin line), T = 170 MeV: S Q = S H (dashed 
line), 7 9 = "/ q r (thin dashed line) and y q = 1 (solid line with 
dot marker); for T = 140 MeV, S Q = S H (dash-dot line) 



Equilibrium yield is much larger than non-equilibrium for 
T < 230 MeV when s/S = 0.04 and S H = S Q , and for 
T < 190 MeV when s/S = 0.03 and S H = S Q . In case 
7 g = 7q r , the yield of ccq is always smaller than equilib- 
rium. The yield of ccs baryons has similar suppression, 
but it becomes larger than equilibrium for smaller tem- 
peratures and yield enhancement for higher T is larger for 
S H = S® then in case of ccq because of large number of 
strange quarks. 

In the figure l23l we show ratios ccq/ J/9 (upper panel) 
and ccs/J/9 (lower panel) as a function of hadroniza- 
tion temperature. These ratios do not depend on dV/dy. 
ccq/ J/9 oc 7 9 does not depend on s/S. For ccq/ J/9 ratio 
we show three cases: chemical equilibrium j s — j q = 1 
(dashed line), S H — S® (solid line) and 7 9 = 7^ r (dash- 
dot line). For ccs /J/9 (ccq/ J/9 oc j s ) we show chemical 
equilibrium case (dashed line), s/S — 0.04: S H = S® 
(solid line with point marker) and j q — j q r (thin dash-dot 
line); s/S = 0.03: S H = S Q (solid line) and lq = 7^ (thin 
dash-dot line). The overall all yields of double charmed 
(strange and non-strange) baryons and anti-baryons is 
clearly larger than the yield of J/9. 



7 Conclusions 

We have considered here in some detail the abundances of 
heavy flavor hadrons within the statistical hadronization 
model. While we compare the yields to the expectations 
based on chemical equilibrium yields of light and strange 
quark pairs, we present results based on the hypothesis 
that the QGP entropy and QGP flavor yields determine 
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Fig. 21. (Color on line) Be mesons yields as function of T for 
chemical equilibrium case with dV/dy = 600 /m -3 for T = 
200 MeV (the upper panel, dashed line), for s/S = 0.03 with 
dV/dy = 600 /m~ 3 for T = 200 MeV (the upper panel, solid 
line), for chemical equilibrium case with dV/dy = 800 /m~ 3 
for T = 200 MeV (the lower panel, dashed line) for s/S = 0.04 
dV/dy = 800 /m" 3 for T = 200 MeV (lower panel, solid line) 



the values of phase space occupancy 7? i = q, s, c,b, which 
are of direct interest in study of the heavy hadron yields. 

For highest energy heavy ion collisions the range of 
values discussed in literature is 1 < 7J 1 < 1.65 and 0.7 < 
7^/7^ < 1.5. However 7^ and 7^ values which are much 
larger than unity arise. This is due to the need to describe 
the large primary parton based production, and consider- 
ing that the chemical equilibrium yields are suppressed by 
the factor exp(— m/T). 

Our work is based on the grand canonical treatment 
of phase space. This approach is valid for charm hadron 
production at LHC, since the canonical corrections, as we 
have discussed, are not material. On the other hand, even 
at LHC the much smaller yields of bottom heavy hadrons 
are subject to canonical suppression. The value of the pa- 
rameter 7^ obtained at a fixed bottom yield iV&, using 
either the canonical, or the grand canonical methods, are 
different, see e.g. Eq. (15) in [33 . Namely, to obtain a given 
yield Nf, in canonical approach, a greater value of 7^ is 
needed in order to compensate the canonical suppression 
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Fig. 22. (Color on line) ccq/N^ (upper lines in each panel) and 
ccs/Nc (lower lines in each panel) baryon yields as a function 
of T. Upper panel: chemical equilibrium case with dV/dy = 800 
fm -3 for T = 200 MeV(dashed line), s/S = 0.04 with dV/dy = 
800 /m~ 3 for T = 200 MeV: S H = S Q (solid line) and 7, = 
7q r (dash-dot line); and lower panel: chemical equilibrium case 
with dV/dy = 600 fm~ 3 for T = 200 MeV (dashed line), and 
s/S = 0.03 S Q = S H (solid line) and 7 9 = 7^ (dash-dot line). 



effect. However, for any individual single-6 hadron, the rel- 
ative yields, i.g. B/B s do not depend on 7^ and thus such 
ratios are not influenced by canonical suppression. More- 
over, as long as the yield of single-6 hadrons dominates the 
total bottom yield: Nb — B + B s + A\, + . . ., also the Nb 
scaled yields of hadrons comprising one 6-quark i.e. ratios 
such as B/Nb, B s /Nb, B c /Nb, etc, are not sensitive to the 
value of 7^ and can be obtained within either the canon- 
ical, or grand canonical method. On the other hand for bb 
mesons and multi-6 baryons the canonical effects should 
be considered. Study of the yields of these particles is thus 
postponed. 

We address here in particular how the yields of heavy 
hadrons are influenced by 7^ /7J 1 7^ 1 and j q 7^ 1. The ac- 
tual values of 7^/7^ we use are related to the strangeness 
per entropy yield s/S established in the QGP phase. Be- 
cause the final value s/S is established well before hadroniza- 
tion, and the properties of the hadron phase space are well 
understood, the resulting 7^ / 7J 1 are well defined and turn 
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Fig. 23. (Color on line) ccq/. J/9 (upper panel) and ccs/J/I' 
(lower panel) ratios as a function of T. Upper panel: chemical 
equilibrium case (dashed line), S H = (solid line) and 7, = 
7q r (dashed-dot line); and lower panel: chemical equilibrium 



case with (dashed line), 
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out to be quite different from unity in the range of temper- 
atures in which we expect particle freeze-out to occur. We 
consider in some detail the effect of QGP hadronization 
on the values of 7^ and 7J 1 . 

One of first results we present (figure [2]) allows a test 
of the statistical hadronization model for heavy flavor: we 
show that the yield ratio cc ss/ (cs cs) is nearly indepen- 
dent of temperature and it is also nearly constant when 
the 4> is allowed to freeze-out later (figure [3]), provided 
that the condition of production is at the same value of 
strangeness per entropy s/S. 

We studied in depth how the (relative) yields of strange 
and non-strange charmed mesons vary with strangeness 
content. For a chemically equilibrated QGP source, there 
is considerable shift of the yield from non-strange D to the 
strange D s for s/S — 0.04 expected at LHC. The expected 
fractional yield D s /N c ~ B s /Nb ~ 0.2 when one assumes 
7? — Tg 1 — 1) the expected enhancement of the strange 
heavy mesons is at the level of 30% when s/S — 0.04, and 
greater when greater strangeness yield is available. 

A consequence of this result is that we find a rela- 
tive suppression of the multi-heavy hadrons, except when 
they contain strangeness. The somewhat ironic situation 
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is that while higher beyond chemical equilibrium charm 
QGP yield enhances production of cc states, beyond chem- 
ical equilibrium enhanced light quarks and strangeness 
multiplicity suppresses this almost by that much. This 
new phenomenon adds to the complexity of interpreta- 
tion of hidden charm meson yield. On the other hand, the 
yield oicc/N^ ~ 210 -3 is found to be almost independent 
on hadronization temperature in case which conserves en- 
tropy at hadronization. We don't know exactly equation 
of state in QGP and so the value of 7 9 which is needed to 
conserve the entropy may be different. If 7 g is larger for 
higher temperatures, suppression of cc is larger for a fixed 
s/S. The same result is found for B c « 5 - 6 l(T 4 iV c 7V h , 
that yield remains considerably larger (by a factor 10 — 
100) compared to the scaled yield in single nucleon mi- 
cleon collisions. 

We have shown that the study of heavy flavor hadrons 
will provide important information about the nature and 
properties of the QGP hadronization. The yield of Bc(6c) 
mesons remains enhanced while the hidden charm cc states 
encounter another suppression mechanism, compensating 
for the greatly enhanced production due to large charm 
yield at LHC. 
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